Interpreting the Wide Scattering of Synchronized Traffic Data by Time Gap Statistics 
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Based on the statistical evaluation of experimental single-vehicle data, we propose a quantitative 
interpretation of the erratic scattering of flow-density data in synchronized traffic flows. A correlation 
analysis suggests that the dynamical flow-density data are well compatible with the so-called jam line 
characterizing fully developed traffic jams, if one takes into account the variation of their propagation 
speed due to the large variation of the netto time gaps (the inhomogeneity of traffic flow). The form 
of the time gap distribution depends not only on the density, but also on the measurement cross 
section: The most probable netto time gap in congested traffic flow upstream of a bottleneck is 
significantly increased compared to uncongested freeway sections. Moreover, we identify different 
power-law scaling laws for the relative variance of netto time gaps as a function of the sampling size. 
While the exponent is —1 in free traffic corresponding to statistically independent time gaps, the 
exponent is about —2/3 in congested traffic flow because of correlations between queued vehicles. 



Recently, the statistical physics and non-linear dynam- 
ics of driven many-particle systems have been key disci- 
plines in the discovery, interpretation, and simulation of 
phenomena in traffic flows [1—3]. The observed instability 
mechanisms, jamming, segregation, breakdown, and clus- 
tering phenomena are now viewed as a paradigm for sim- 
ilar phenomena in granular and colloidal physics [3-5] , in 
biology (social insects) [6] , logistics (instability of supply 
chains), and economics (business cycles) [7,4]. Empirical 
investigations [8-11] have particularly stimulated the de- 
velopment of traffic theory [11-13]. For example, it was 
possible to identify constants of traffic dynamics such as 
the propagation speed C w — 15 km/h of wide jams [14]. 
According to analytical calculations by Kerner et al. for 
a macroscopic traffic model assuming uniformly behav- 
ing driver- vehicle units with identical parameters, fully 
developed wide traffic jams should be characterized by a 
self-organized flow-density relation 
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where p denotes the vehicle density, while the average 
netto time gap T (the time clearance) and the maxi- 
mum density p max are assumed to be fixed parameters 
[15]. The dependence of the flow J on the density p in 
Eq. (1) is sometimes referred to as "jam line" J(p) and 
does not necessarily agree with the high-density part of 
the steady-state flow-density relation Qi(p) for station- 
ary and homogeneous traffic flow, the so-called funda- 
mental diagram [15]. We note, however, that a linear 
flow-density relation for wide moving jams is not yet 
fully confirmed by empirical measurements and that ag- 
gregation methods used to determine macroscopic traffic 
data from single- vehicle data have a non-trivial impact on 
the measurement results [16,17]. Nevertheless, the slope 
C = dJ(p)/dp — — l/(pmaxT 1 ) is usually identified with 
the average propagation velocity of wide moving jams 



[15]. In contrast to wide moving jams, synchronized traf- 
fic flow is a form of congested flow which is mostly found 
upstream of inhomogeneities (e.g. freeway bottlenecks) 
and claimed to show a completely different behavior [9]. 
It is characterized by an erratic motion of time-dependent 
flow-density data in a two-dimensional area (and a syn- 
chronization of the time-dependent average vehicle veloc- 
ities among neighboring lanes) [8]. More specifically, in 
synchronized flow, the slopes 

C(t k+1 ) = {Q(t k+1 ) - Q(t k )]/[p(t k+1 ) - P (t k )} (2) 

of the lines connecting flow-density data measured at a 
given freeway cross section at successive times t k and t k +\ 
erratically take on positive and negative values, charac- 
terizing a complex spatio-temporal dynamics [8]. This is, 
in fact, one of the most controversial subjects in traffic 
theory. It has not only been the reason for Kerner's fun- 
damental criticism of all traffic models assuming a fun- 
damental (steady-state) relation Qf(p) between the flow 
Q and the density p (i.e. of the vast majority of models 
suggested up to now) [18,1]. It has also triggered a flood 
of publications in physics journals with various sugges- 
tions how to describe this wide scattering. The proposed 
interpretations reach from shock waves propagating for- 
ward or backward with speed C(t) [8,18], changes in the 
driving behavior in response to the traffic situation (in- 
cluding "frustration effects" ) [19,17], anticipation effects 
[20,21], non-unique solutions [22], a trapping of vehicles 
[11], or multiple metastable oscillating states [13]. Non- 
unique solutions are also expected for car-following mod- 
els, in which the vehicle acceleration dv/dt is not a unique 
function of the distance or velocity [23] . Other explana- 
tions based on a mixture of different vehicle types such as 
cars and trucks [24] or a heterogeneity of time headways 
[25] are easier accessible to empirical verification. In the 
following, we will therefore focus on these. The particu- 
lar approach of this study is its restriction to empirical 
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analyses to avoid theory-driven interpretations. Although 
it was been seriously questioned that relation (1) would 
also be applicable to flow-density data of synchronized 
congested traffic flow, by correlation analysis we will sur- 
prisingly find a strong empirical support for the validity 
of relation (1), if the average netto time gap is treated as 
a dynamical variable. 

This requires single-vehicle data, which were recorded 
by double induction-loop detectors for the Dutch freeway 
A9. The measured data include, for each lane, the pas- 
sage time of vehicle i, its velocity Vi, and its length k. 
From this, we determine the individual netto time gaps 
as Tj = t° — t°_ 1 — k/vi. As discussed in Ref. [27], we 
obtain macroscopic quantities by averaging over N = 50 
successive vehicles, which gives a better statistics than 
averaging over a fixed time interval. Thus, we obtain the 
empirical traffic flow at time t k 



Q(t k ) = i/(t°-tl 1 ) = N/(tl N 

age netto time gap by 
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Likewise, we define the vehicle density p and the maxi- 
mum density p max by 
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For a given loop detector on the left lane, we con- 
sider all time intervals with congested traffic (defined by 
p(tk) > 45 vehicles per kilometer and lane) and com- 
pare the temporal changes AQ(' £+1 ' = [Q(t k +i) — Q{tk)] 
of the empirically measured flow with the changes pre- 
dicted by Eq. (1). We will compare three hypotheses by 
correlation analysis: (i) p is treated as an independent 
variable defined by (4), while T and p max are treated 
as parameters, (ii) p and 1/T are independent variable 
defined directly from the single-vehicle data via (4) and 
(3), while p max is a parameter, (iii) p, 1/T, and l/p max 
are independent variables defined by (3) to (5). The first 
order Taylor approximations for the temporal change 
AJ (a) F1) = [ J (a)(tk+i) - J( a )(tk)} of the corresponding 
theoretical vehicle flows are obtained by differential of 
relation (1) with respect to the independent variables, 
but not the parameters: 

[p(t k+1 ] - p(t k )} , (6) 
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Our statistical analysis gives the following correlations 
Corr(AQ,AJ (Q) ) = (£ fc AQ« AJ ( ( Q fe j)/[(£ fc AQ« 2 ) 

(Et ^-^(a) 2 )] 1 ' 2 f° r tnc hypotheses a =i, ii, and iii: 
Corr(AQ, AJ (i) ) = 0.347, Corr(AQ, A J (ii) ) = 0.918, 
Corr(AQ, A J(iii)) = 0.938. As a consequence, hypothesis 
(i) assuming a linear variation of the flow with the density 
is in fact a poor description of "synchronized" congested 
traffic data due to their two-dimensional scattering. How- 
ever, hypothesis (ii) taking into account the dynamical 
variation of the average netto time gap T yields a strong 
correlation, which indicates that we have identified the 
main reason for the wide scattering (see Fig. 1). Taking 
into account a variation of p max with the truck fraction 
improves the correlation coefficient only a little. Hence, 
the maximum density, which influences the density-offset 
of the jam line, is an unimportant explanatory variable. 
In contrast, the netto time gap T determining the slope 
of the jam line (see Eq. (1)) is highly significant. For this 
reason, we will investigate the characteristic features of 
the netto time gap distribution more closely. 

The variation of the average time gaps with time can 
only be relevant, when it is large. The time gap distri- 
bution is, in fact, surprisingly wide, and it has so-called 
heavy tails (see Fig. 2). To investiagte this, let us calcu- 
late the variance of the sample-averaged time gaps as a 
function of the sampling size N, 
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where (j + 1) runs over the first vehicle indices of all 
possible samples of N consecutive vehicles in the data 
set of size I, = X^i=f+i ^ s defined as moving 

average of the time gaps of the N next vehicles, and T = 
J2l = iTi/I = Tq^ is the global average. To account for 
artifacts caused by the daily variation of traffic flow, we 
applied a high-pass filter (with cutoff period N c = 500) to 
the single- vehicle data before calculating Eq. (7), which 
limits the scaling range of relation (8). For a given N, the 
variance (7) decreases with increasing density (see Figs. 2 
and 3), as less and less space is available for time gaps 
larger than the safe time gap, but the variance is also a 
function of the sampling size N . For free traffic flow (with 
p < 15 veh/km/lane), we observe the power law behavior 
Var(T) oc 1/N, as expected for statistically independent 
time gaps in free traffic. At high vehicle densities (with 
p > 35 veh/km/lane), however, we find the power law 



Var(T) oc iV 7 



(8) 



with an exponent 7 w —0.67 (see Fig. 3). That is, the rel- 
ative variance decreases much slower with the sampling 
size than expected, implying that the time gaps do not 
average well and fluctuations of average time gaps re- 
main significant for reasonable sampling sizes N. This is 
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related with surprisingly wide time gap distributions and 
results from correlations between queued vehicles, which 
are probably due to platoon formation [11], but may also 
be caused by dynamical vehicle interactions. The repro- 
duction of this scaling law will be a serious test for mi- 
croscopic traffic simulation models, which are needed to 
reveal the detailed mechanism behind it. 

According to Figure 2, the time gap distribution is not 
only dependent on the density, but also on the measure- 
ment cross section (more specifically: on the relative lo- 
cation with respect to bottlenecks of the freeway). It has 
been suggested [26] that the increase of the most prob- 
able brutto time gap (time headway) At? = t® — t^_ 1 — 
Ti + U/vi in congested traffic compared to free traffic 
[10,27] is due to the drop of the vehicle velocity Vi, in- 
creasing h/vi. However, an increase is also observed for 
the most probable netto time gaps Tj (see Fig. 2). This 
increase is most pronounced at cross sections upstream 
of bottlenecks (e.g. on-ramps), where vehicles are queuc- 
ing. We interpret this increase of the netto time gaps as 
"frustration effect" of drivers after a considerable queu- 
ing time, which is supported by other observations as 
well [3,28]. The frustration effect is most significant at 
freeway sections, where the accumulated time of driving 
under congested conditions is high, and it decays with 
increasing recovery time. At freeway sections, where seri- 
ous congestion never occurs, there is no relevant increase 
in the most probable netto time gap with growing density 
(see Fig. 2). Finally, we note that the time gaps in front 
of long vehicles ("trucks") are much higher, on average, 
than in front of short vehicles ( "cars" ) , as expected (see 
Fig. 2). This supports the theory suggested in Ref. [24]. 

In conclusion, we have found an interpretation of the 
wide scattering of flow-density data in "synchronized" 
congested traffic based on the jam relation J(p), but tak- 
ing into account the time-dependent variation of empir- 
ical netto time gaps among successive cars. This varia- 
tion is related to time-dependent changes of the slope 
C = — l/(pmax? 1 ) of the jam line and causes the two- 
dimensional scattering together with the variation of the 
vehicle density. The surprisingly strong variation of the 
average time gaps was due to the fact that the power law 
relating the variance of the time gaps with the sampling 
size had a considerably smaller exponent (7 w —2/3) 
than expected (7 = —1). Hence, distinguishing differ- 
ent forms of congested traffic based on the scattering is 
questionable. We also found that the increase in the most 
probable time gap in congested traffic is site-dependent 
and not attributed to a drop in the vehicle velocity, indi- 
cating a "frustration effect". Traffic models considering 
the dynamical changes of the netto time gaps can repro- 
duce empirical observations more realistically [24,17]. 
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FIG. 1. (Color online)The two-dimensional scattering of 
empirical flow-density data in synchronized traffic flow of 
high density p > 45 veh/km/lane (see (a)) is well repro- 
duced by the jam relation (1), when not only the variation 
of the density p, but also the empirically measured varia- 
tion of the average time gap T and the maximum density 
Pmax is taken into account (see (b)). The similarity between 
the experimental data and relation (1) is partly because the 
density p(t&) plotted in (a) and (b) (the x-value) is de- 
termined with the same formula (4), but the agreement of 
the empirical flow Q — l/(Ti + h/t>i) and of the theoret- 
ical relation J(p,T,p max ) = [1 - (Z;)/<Vi(T; + k/vi))}/{Ti) 
(the y-values) is not trivial, as even low-order approxima- 
tions of these formulas differ. The pure density-dependence 
J(p) (thick black line) is linear and cannot explain a 
two-dimensional scattering. However, variations of the aver- 
age time gap T change its slope — l/(pmaxT) (see arrows), 
and about 95% of the data are located between the thin lines 
J{p,T± 2AT, 1/7) = (1 - pl)/(T± 2AT), where I = 3.6m is 
the average vehicle length, T = 2.25 s the average time gap, 
and AT = 0.29 s the standard deviation of T. The predicted 
form of this area is club-shaped, as demanded by Kerner [18]. 
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FIG. 2. (Color online) (Netto) time gap distributions at 
different cross sections of the Dutch freeway A9 (see the 
sketch in subfigure (i)), separately measured for the fast and 
slow lane, for car-car, car-truck, and truck-car interactions, 
and for different density regimes, namely free flow (p < 10 
veh/km/lane), low densities (free traffic with 10 < p < 25 
veh/km/lane), medium densities (mostly congested traffic 
with 25 < p < 50 veh/km/lane), and high densities (con- 
gested traffic with p > 50 veh/km/lane, only observed up- 
stream of serious bottlenecks). Densities were determined as 
50- vehicle averages, and trucks were defined by a vehicle 
length U > 6 m. The time gaps of trucks with respect to cars 
are, on average, considerably higher than the ones of cars with 
respect to other vehicles, which causes a large variation of in- 
dividual time gaps. However, even the variation of time gaps 
among cars is considerable. A remarkable point is the increase 
in the most probable time gap from 0.75 s in free traffic (on 
the fast lane) to considerably higher values in congested traf- 
fic upstream of a bottleneck (1.2 s at medium densitities and 
1.9 s at high densities at Loop 4 compared to 0.75 s at Loop 12 
in all density regimes and 1 s at Loop 7 at medium densities, 
see (a), (c), (e); the values for the slow lane are somewhat 
higher). This "frustration effect" is still slightly active imme- 
diately downstream of a bottleneck (Loop 7), in contrast to 
freeway sections, which are never seriously congested (Loop 
12). It is also observed for finite minimum bumper-to-bumper 
distances of, for example, lo = 1.5 or 3 meters, which is con- 
sidered by replacing U by (h + lo) everywhere. 
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FIG. 3. (Color online) Variance of the average time gaps T 
as a function of the sampling size N, separately measured for 
free and congested flows at loop 4 on the fast lane. While in 
free flow, the relative variance is consistent with the expected 
power law exponent of 7 = —1, it decreases much slower with 
N in the congested traffic regime, namely with an exponent 
of 7 ss —0.67, causing the significant variation of T. 
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